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Abstract. Conventional astrophysical observations have failed to provide 
stringent constraints on physical processes operating in the interior of the stars. 
However, satellite missions now promise a solution to these problems by provid- 
ing long-term high-quality continuous data which will allow the application of 
seismic techniques. With this in mind, and using the Sun as our astrophysical 
laboratory, our aim is to determine if Corot- and Kepler-like asteroseismic data 
can constrain physical processes like microscopic diffusion. We test to what ex- 
tent can the observed atmospheric abundances coupled with p-mode frequencies 
safely distinguish between stellar initial chemical composition and diffusion of 
these elements. We present some preliminary results of our analysis. 



1. Introduction & Method 

If we have the following observations for a solar-type star: a set of p-mode 
frequencies, measured abundances, and an identified g mode, we can investigate 
how these observations may constrain some of the physical processes occuring 
in the interior of a star. Suppose that the stellar model is a main sequence 1 
Mq star, and that we can calibrate the age of the star using the radius and 
luminosity. Then the free parameters that remain to be fitted in the star are 
the initial mass fractions of hydrogen and metals Xq and Zq, and the velocities 
of diffusion of these elements Vx, Vy and Vz, (assuming that other microscopic 
processes may not be "observable"), where Y denotes helium. Let's denote these 
set of parameters as P. Using P we can calculate a stellar model which will 
give us the expected observables, B (expected p-mode frequencies, abundances 
etc.). If we have a set of observations O, then in order to find the P that most 
adequately reproduce O, we can use the following equation iteratively, until we 
make O as close as possible to B: 

5P = VW-iU^^^. (1) 
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Figure 1. The relative uncertainties as a fraction of the initial values of 3%, 
24%, 189%, 217%, 422% (details in text) as extra observables are included. 



SP are the parameter corrections to make to the initial guess Pq, e are the 
measurement errors, and UWV^ is the singular value decomposition of the 



sensitivity matrix, which can be written as / e a nd W ^ has some zero values 



to stabilize the inversion fsee ICreevevI (120081 . 120091 ')) . 

Once P have been found, then the theoretical uncertainties <t(P) can be 
calculated for each j = 1, 2, ...A^ parameter 
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The partial derivatives were calculated using the CESAM c ode Code d'E volutio n 
SteUaire Adaptatif et Modulaire f.MorelllQQ? ) and ADIPLS ( Christensen-DalsgaardI 
20071') ■ We took the stari dard solar model as a reference using the abundances of 
Grevesse &: Noelsl ( 19931 ). and calibrated t hese models iii radiu s and luminosity 
to an accuracy of ~ 10~^ as described in iMathur et all (l2007l ). The reference 
parameters for the models are {Xq,Zq) = (0.708,0.0195) and {Vx,Vy,Vz) = 
(2.307, -7.063, -5.635) x 10"^ cm s'K 



Results and Discussion 



2.1. Impact of each observable on the uncertainties 

Using only p-mode frequencies to determine the parameter uncertainties (using 
Eq. [2]) the diffusion coefficients Vx, Vy and Vz can not be disentangled. Their 
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uncertainties remain too large (189%, 217% and 422% respectively) when we 
assume that we have only these observations with errors of the order of 1.3/iHz, 
meanwhile the values for (t(Xo) and (t{Zq) are 3% and 24%. However, to show 
the effect that each of the observations has on the determination of the pa- 
rameters, Figure 1 shows how the uncertainties in each of the parameters reduce 
relative to these values, as extra observables such as abundances and g modes are 
included. Each of the bars represents from left to right the following additions 
to the set of observations: 

• Blue reduce the errors on the p-mode frequencies by almost a factor of 2 
(from 1.3 /iHz - 0.7/iHz) 

• Green include abundance measurements of He, C, N and O 

• Yellow include abundance measurements of Li and Be (e = 0.1 dex for all 
abundances) 

• Red include 1 identified g mode 

After including all of these observations, the final parameter uncertainties are: 
1%, 6%, 28%, 57% and 72% respectively. 

It is interesting to note the various effects of adding in new observables. 
Firstly, reducing all of the uncertainties by almost a factor of 2 should result 
in a corresponding reduction in the uncertainties by this same amount. The 
blue bars clearly indicate that this is the case. When we include the abundance 
measurements of He, C, N, O with an observational error of 0.1 dex, the uncer- 
tainties in both Xq and Zq reduce by almost a factor of two (green). These same 
abundances have little effect for constraining the diffusion velocities. However, 
including Be and Li measurements has a significant impact on the determina- 
tion of the metal diffusion velocity Vz (yellow), reducing its uncertainty to 72% 
(~20% of its original value). Finally, the inclusion of one identified g mode has 
most impact on constraining Vx and then Vy (red). 

2.2. Recovery of original parameters 

The final uncertainties may be sufficiently low that the observations do contain 
enough information to be able to distinguish safely between initial metal mass 
fraction Zq and diffusion of this Vz (our original scientific question). Assume 
that we have two new sets of observations, and that these O really come from a 
solar model with 1) lower Zq and 2) lower Vz- If we use Eq. [T]to fit the new O 
(separately) while using the reference set of parameters as Pq (the initial guess), 
we will obtain non-zero negative values of 6Zq and 5Vz respectively, only if the 
observations contain sufficient information. If there is not enough information 
contained in the observations, then 5Zq ~ 5Vz ~ 0. 

Figure [2] shows the values of 5P (relative to the original value of P, for 
scaling purposes) for some of the parameters of the models, when attempting to 
fit a set of observations generated from a model with a reduced Zq (red) and a 
reduced Vz (blue). As the figure shows, some non-zero 6P are needed to account 
for the new observations. 

• When inverting using the lower Zq model observations (red) , 5Zq is nega- 
tive. This means that we need to reduce the original Zq value to a lower 



4 



Inversion from reference model 



100 
90 
80 
% 70 



60 

50 
40 



• -Vj model 

• -Z model 



Figure 2. The parameter changes SP needed to make to Pq to correctly fit 
the new observations. These new observations were generated using a model 
with a lower Zq (red) and a lower Vz (blue). 



one — consistent with the input model. However, we also see a nega- 
tive contribution to Vz; ideally Vz would not change, because it was not 
changed in this model, just as Xq does not change. 

• When inverting using the lower Vz model observations (blue), we obtain 
a much more negative 5Vz (but no change in Zq), indicating a necessary 
downward revision of only Vz — consistent with the input model. 

These results are preliminary, but encouraging: there is some indication in 
the set of observables that we may be able to distinguish between initial metal 
abundance and diffusion of this using abundance measurements and seismic data. 
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